Titanium-decorated carbon nanotubes: a potential high-capacity hydrogen storage 
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We report a first-principles study, which demonstrates that a single Ti atom coated on a single- 
walled nanotube (SWNT) binds up to four hydrogen molecules. The first H2 adsorption is dissocia- 
tive with no energy barrier while other three adsorptions are molecular with significantly elongated 
H-H bonds. At high Ti coverage we show that a SWNT can strongly adsorb up to 8-wt% hydro- 
gen. These results advance our fundamental understanding of dissociative adsorption of hydrogen 
in nanostructures and suggest new routes to better storage and catalyst materials. 
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Developing safe, cost-effective, and practical means of 
storing hydrogen is crucial for the advancement of hy- 
drogen and fuel-cell technologies 0. The current state- 
of-the-art is at an impasse in providing any material that 
meets a storage capacity of 6-wt% or more required for 
practical applications [l[ 0, 0, H, @, 0- Here we 
report a first-principles computation of the interaction 
between hydrogen molecules and transition metal atoms 
adsorbed on carbon nanotubes. Our results are quite 
remarkable and unanticipated. We found that a single 
Ti-atom adsorbed on a SWNT can strongly bind up to 
four hydrogen molecules. Such an unusual and complex 
bonding is generated by the concerted interaction among 
H, Ti, and SWNT. Remarkably, this adsorption occurs 
with no energy barrier. At large Ti coverage we show 
that a (8,0) SWNT can store hydrogen molecules up to 
8-wt%, exceeding the minimum requirement of 6-wt% for 
practical applications. Finally, we present high temper- 
ature quantum molecular dynamics simulations showing 
that these systems are stable and indeed exhibit associa- 
tive desorption of H2 upon heating, another requirement 
for reversible storage. 

Recent experiments [il Ho| and calculations [TfL H2I H3I 
suggest that it is possible to coat carbon nan- 
otubes unifor mly with Ti atoms without metal segrega- 
tion problems [13. Here we show that such Ti-coated car- 
bon nanotubes exhibit remarkable hydrogen absorption 
properties. Below we will present our results in detail 
for a (8,0) nanotube and briefly for four armchair (n,n) 
(n=4,5,6, and 7) and five zigzag (n,0) (n=7,8,9,f0,lf , and 
12) nanotubes. 

A single Ti atom on an (8,0) SWNT has a magnetic 
ground state with S=f and a binding energy of 2.2 eV; 
this will serve as our reference system, denoted t80Ti. 
In order to determine different reaction paths and prod- 
ucts between H2 and t80Ti, we have carried out a series 
of single-energy calculations as H2 molecules approaches 
t80Ti and when there are large enough forces acting 
on H2 molecules, we let the atoms evolve according to 
the quantum mechanical forces obtained from density 



functional theory (DFT) calculations 16]. We used the 
conjugated-gradient (CG) minimization and optimized 
both the atomic positions and the c-axis of the tube. 

The energy calculations were performed within the 
plane-wave implementation of the generalized gradi- 
ent approximation (GGA-PBE)H3 to DFT. We used 
Vanderbilt ultra-soft pseudopotentials^^ treating the 
following electronic states as valence: Ti: 3s 2 3p 6 3d 2 4s 2 , 
C: 2s 2 2p 2 and H: Is. The Monkhorst-Pack special k- 
point scheme^] is used with 0.025 A -1 k-point spacing 
resulting in 5 k-points along the tube axis. The cutoff 
energy of 350 eV is found to be enough for total en- 
ergies to converge within 0.5 meV/atom. The calcula- 
tions are carried out in a tetragonal supercell geometry 
of 20A x 20A x c where c is taken to be twice the lattice 
constant of SWNT along its axis. 

Figure la shows the energy variation from non-spin- 
polarized calculations as a single H2 molecule approaches 
t80Ti. The energy first decreases slowly as the hydro- 
gen gets closer to the nanotube and Ti. However, as the 
charge overlap gets large, the H 2 molecule is attracted 
towards the Ti atom with a sudden decrease in the en- 
ergy. At this point, the H2 molecule is still intact with 
a significantly increased H-H bond length of 0.86 A. The 
second sudden decrease in energy is achieved by dissoci- 
ating the H2 molecule into two H atoms. At this point, 
the H-H distance increases from 0.86 A to 2.71 A. The in- 
teraction between H2 and t80Ti is always attractive and 
therefore H2 is absorbed onto a Ti atom without any en- 
ergy barrier. The final geometry is shown in the inset to 
Fig. la, with relevant structural parameters given [20j. In 
order to calculate the binding energy for this dissociative 
adsorption, we calculate the total energies of the initial 
t80Ti and H2 state (dashed line in Fig. la) and the final 
t80TiH2 state (dotted line in Fig. la) from spin-polarized 
calculations. We obtain the binding energy to be 0.83 eV 
(Fig. 1). We note that spin-polarized calculation lowers 
the total energy by about 0.6 eV with respect to non- 
spin-polarized calculations and yields a triplet magnetic 
ground state (i.e. S=l) for the initial t80Ti and H2 sys- 
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FIG. 1: Energy versus reaction paths for successive dis- 
sociative and molecular adsorption of Eh over a single Ti- 
coated (8,0) nanotube. {a)H 2 + t80Ti -> t80TiH 2 . (b) 
2H 2 +t80TiH 2 -> t80TiH 2 -2H 2 . (c) H 2 + t80TiH 2 -2H 2 -> 
t80TiH 2 — 3H 2 . The zero of energy is taken as the sum of the 
energies of two reactants. The relevant bond distances and 
binding energies (i?s) are also given. 



torn. However once the hydrogen molecule is attached to 
Ti, the system is non-magnetic and spin-polarized calcu- 
lations are not necessary. 

Remarkably, it is also energetically favorable for 
the TiH2 group to complex with additional hydrogen 
molecules. Fig. lb shows the energy variation as two hy- 
drogen molecules approach the Ti atom; one from each 
side of the TiH2 group. As in the case of single adsorp- 
tion, the energy always decreases, first slowly and later 



very rapidly at which point both hydrogen molecules are 
strongly attached to the t80TiH2 system. We denote 
the final product as t80TiH2-2H2, which is shown in 
Fig. lb. The total energy change upon adsorption is 
about 0.89 eV (i.e. 0.45 eV/H 2 ). Unlike the first adsorp- 
tion, the two hydrogen molecules are in intact but with 
a rather elongated bond length of 0.81 A. This 10% in- 
crease is rather reminiscent of the elongated H-H bonds 
observed in the metal-dihydrogen complexes first synthe- 
sized by KubasjUllil. 

Fig. lc shows the energy evolution when a fourth hy- 
drogen molecule approaches the t80TiH2-2H2 system 
from the top. The energy again decreases continuously, 
indicating a zero-energy barrier. The final product, de- 
noted as t80TiH 2 -3H 2 , is shown in the inset with the 
relevant structural parameters. The energy gained by the 
fourth adsorption, 0.34 eV/H 2 , is slightly smaller than 
for the other cases but is still substantial. The H-H dis- 
tance of the top H2 is 0.86 A. Several attempts to add 
a fifth hydrogen molecule at a variety of positions failed, 
suggesting a limit of 4 H2/TL 

The final optimized structures shown in Fig. 1 need not 
be the global minimum. Among many different isomers 
tried for the four-H2 system, we found a very symmet- 
ric configuration (denoted as t80Ti-4H 2 ) that is 0.1 eV 
lower in energy than t80TiH 2 -3H 2 (Fig. 2a). Here all 
four hydrogen molecules stay intact and benefit equally 
from bonding with the Ti atom. The average H-H bond 
distance is about 0.84 A and each molecule has an ex- 
cess charge of about 0.15 e. The average binding en- 
ergy per H2 molecule is 0.54 eV, which suggests that the 
bonding is an unusual combination of chemisorption and 
physisorption. 

To test system stability, we performed high- 
temperature quantum molecular dynamics (MD) simu- 
lations on t80TiH2-3H2 [13 • I n total, we have carried 
out a 1.5 ps MD simulations with a Langevin thermo- 
stat for temperatures ranging from 200 K to 900 K. We 
observe an associative desorption of a H2 molecule (i.e, 
two H atoms come together to form H2 and leave the 
system) at around 800 K^. While a 1.5 ps time MD 
simulation is already computationally quite costly, it is 
not long enough to get statistically meaningful values for 
the temperature. However, it does suggest that the sys- 
tem is quite stable and that it is possible to extract the 
H2 molecules without breaking the Ti-C bonding or re- 
moving the T1H2 from the nanotube. 

Projecting the plane waves on the pseudo-atomic or- 
bitals (Fig. 2b) gives detailed information on the nature 
of the bonding. The density of states indicates that Ti 
3s— and 3p-orbitals are mostly intact and lie below - 
20 eV. In the energy range of from -10 eV to zero, Ti 
contributes only d-electrons while Ti 4s-electrons are al- 
most absent, probably promoted to Ti d-orbitals. Fig. 2b 
shows that the binding state just below Ep has a major 
contribution from Ti d-orbitals, together with the carbon 
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FIG. 2: (a) Two different views of the optimized structure 
of t80Ti-4H2. The relevant structural parameters are: H- 
H=0.84 A, Ti-H=1.9 A, Ti-C=2.17 A, Ti-C'=2.4 A. (b) The 
projected density of states (PDOS) at T-point contributed 
from Ti, four H2 molecules and the six carbons of the hexagon 
on which Ti and H2 molecules are bonded, (c) the a*- 
antibonding orbital of four H2 complex, (d-f) isosurface of the 
state just below Ef at three different values; at ^ = 0.08 it 
is mainly Ti-d orbital, at >]/ = 0.04 the hybridization between 
d-orbital, two carbon pi-orbital and 4H2 cr*-antibonding is 
apparent. At \P = 0.02 it is clear that the other four carbon 
atoms are also involved in the bonding. 



p-orbitals and hydrogen s-orbitals. In Fig. 2, we showed 
the isosurface plot of this state. From this figure it is 
clear that a filled Ti d-orbital is hybridized simultane- 
ously with the er*-antibonding molecular orbital of 4H2 
complex (see Fig. 2c) and the 7r|-antibonding orbital of 
the C6-ring of the SWNT. The Mulliken analysis indicate 
that there is about 1 e transfer to the 4H2 antibonding 
orbital. Similarly PDOS in Fig. 2b shows that the o- 
bonding of the AH2 group forms a band between -10 and 
-6 eV and weakly hybridized with d-orbitals. Integrat- 
ing these peaks over energy and k-points suggests that 
there is about 0.4 e electron transferred to the Ti empty 
d-orbitals. Hence the bonding mechanism for t80Ti-4H2 
seems to be very similar to Dewar, Chatt, and Duncan- 
son modelji^, where the interaction has often viewed as 



a donation of charge from the highest occupied orbital of 
the ligand to the metal empty states and a subsequent 
back donation from filled d-orbitals into the lowest un- 
occupied orbital of the ligand. Finally, the second peak 
around -1 eV in PDOS corresponds to the hybridization 
of one d-orbitals with the tt orbitals of the C6-hexagon of 
SWNT, which is responsible for the bonding of Ti atom 
to the nanotube. 

In summary, our analysis of PDOS and molecular or- 
bitals clearly indicates that we need two occupied d- 
orbitals; one for molecular bonding of the hydrogens and 
the other to bind the metal to the nanotube. We also 
expect that the ionization potential (IP) of the metal 
is important, which controls the amount of back charge 
transfer to the hydrogen antibonding state. When a sin- 
gle H2 molecule is introduced to t80Ti, it seems that Ti is 
able to donate just enough charge to the a* -antibonding 
state, causing dihydrogen to be unstable against disso- 
ciation of H2 (Fig. la). However when more hydrogen 
molecules are added to the system, the charge transfer 
per H 2 molecule is not enough to destabilize the dihydro- 
gen state and therefore the absorption becomes molecular 
(Fig. 2a). 

From above discussion it is clear that we need both 
filled and empty d-orbitals for the metal-hydrogen com- 
plex formation. We tried the same thing with alkali 
and alkaline earth metals, such as Li, Mg and failed. 
Heavy transition metals with diffusive d-orbitals, such 
as Pt and Pd, are also not good candidates for the 
molecular absorption. Such metals are known to inter- 
act with the <7*-antibonding of the hydrogen molecules 
strongly, destabilizing dihydrogen structure against clas- 
sical hydride formation. Our preliminary results based 
on a structural optimization starting from t80TiH 2 -3H 2 
with Ti replaced by Pt/Pd indeed indicate that the side 
hydrogen molecules are not bonded and leave the system 
immediately |3 • We did observe that two H2 do indeed 
bind to Pt/Pd forming a PtHi (or PdH4) classical hy- 
dride cluster, which was not bonded to the nanotube. 
However we expect the light transition metals like Sc 
and V to show similar behavior since they have both oc- 
cupied and empty d-orbitals with a similar IP to Ti. We 
are currently investigating a large number of transition 
metals and the results will be published elsewhere [2^. 

It is important to know if the results reported above 
for (8,0) SWNT hold for other nanotubes and how they 
depends on the chirality and tube radius. Therefore we 
have also studied four (n,n) (n=4,5,6 and 7) and five 
(n,0) (n=7,8,9,10, 11, and 12) nanotubes and details 
will be published elsewhere j24|- Briefly we find that the 
binding energies of TiH 2 , TiH 2 -3H 2 , and Ti-4H 2 groups 
have a weak but complicated dependence on the tube 
radius, conduction and valence band energy levels and 
band gaps. Therefore we did not find a simple trend like 
1/R for the binding energies. This suggests that nan- 
otubes with larger diameter can also show the similar 
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a C 8 TiHg (5.3 wt%) b CJiHg (7.7 wt%) 




FIG. 3: Two high-density hydrogen coverage on Ti-coated 
(8,0) nanotube. 



effect. For the fargest nanotubes that we have studied, 
the binding energy of Ti is reduced to 1.8 and 1.6 eV 
for (12,0) and (7,7), respectively. The binding energy for 
the first H 2 for (7,7) is about 0.66 eV/H 2 , slightly re- 
duced from 0.83 cV/H 2 for (8,0) nanotube. We find that 
those nanotubes with significant band gaps (like (8,0), 
(10,0) and (11,0)) show the strongest hydrogen absorp- 
tion. For example (10,0) has 2.5 eV/4H2 binding energy 
for T1H2-3H2 while the binding energy for (11,0) and 
(8,0) is about 2.1 eV/4H2. We also find interesting small 
differences between (n,n) and (n,0) nanotubes. For ex- 
ample the TiH2-3H 2 configuration (Fig. lc) is the ground 
state for the (n,n) nanotubes, while it is Ti-4H2 (Fig. 2a) 
for (n,0) nanotubes. In conclusion, the phenomena that 
a single Ti atom absorbed on hexagonal phase of SWNT 
can bind up to four molecular hydrogen is a very general 
and novel result and holds for a very large number of nan- 
otubes and possible for other carbon-based nanoclusters 
(i.e. C 6O 0)- 

To this point we have discussed the interaction of 
H2 with a single Ti atom bonded to a nanotube, but 
clearly one can imagine attaching more Ti to a nan- 
otube, thereby increasing the hydrogen storage capac- 
ity. In order to show the feasibility of this approach, we 
present two simple cases where Ti covers j and \ of the 
hexagons. The optimized bond- lengths and other param- 
eters of the structures shown in Fig. 3 are very similar 
to those in the single-Ti case, indicating that the system 
has the capacity to have many Ti and hydrogen. In fact, 
these configurations, which have the chemical formulas, 
CsTiHg and C4T1H8, store approximately 5-wt% and 8- 
wt% hydrogen, respectively. These numbers are based 
on the assumption that t80TiH 2 group will also release 
the hydrogen molecule without difficulty. This is a good 
assumption since the binding energy of H 2 in t80TiH 2 is 
about 1/3 of the binding energy of t80Ti. 

The calculated binding energy^ for CsTiHs is 0.43 
eV per H2, which is somewhat reduced from the average 
binding energy of 0.54 eV in t80Ti-4H 2 . For C 4 TiH 8 , 
the binding energy is further reduced to 0.18 eV. This re- 



duction is due to Ti-Ti interaction which is also respon- 
sible for the increased binding energy of Ti on nanotubes 
(i.e. 2.8 eV/Ti compared to 2.2 eV for a single Ti on a 
(8,0) SWNT). The stability of C 4 TiH 8 is further checked 
by MD simulations on a 1 x 1 x 2 supercell, which did not 
show any Ti metal segregation. Furthermore, we checked 
the stability of the system against releasing one and eight 
of the top hydrogen molecules (see Fig. 3b), both yielded 
a binding energy of 0.13 eV/H2. This is the weakest bond 
in the system and yet it is 4-5 times stronger than the van 
der Walls interactions between hydrogen and SWNt|27|. 

In conclusion, using the state-of-the-art first-principles 
total-energy calculations we have shown that each Ti 
atom adsorbed on a SWNT can bind up to four hydro- 
gen molecules, a remarkable and totally unanticipated 
finding. The mechanism of the bonding is explained 
by a unique hybridization between Ti-d , hydrogen cr*- 
antibonding and SWNT LUMO orbitals (Fig. 2). These 
results advance our fundamental understanding of dis- 
sociative and molecular chemisorption of hydrogen in 
nanostructures, a fundamental step towards novel mate- 
rials needed for hydrogen production, storage, and con- 
sumption in the fuel cells. They also suggest a possi- 
ble method of engineering new nanostructures for high- 
capacity storage and catalyst materials. 
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